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Abstract 

We construct the ?7l(3) x Ur{?>) effective lagrangian which encodes the dy- 
namics of the low energy pseudoscalar excitations in the Color-Flavor-Locking 
superconducting phase of QCD at finite quark density. We include the effects 
of instanton-induced interactions and study the mass pattern of the pseu- 
doscalar mesons. A tentative comparison with the analytical estimate for the 
gap suggests that some of these low energy momentum modes are not stable 
for moderate values of the quark chemical potential /x. 
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I. INTRODUCTION 



The fascinating possibility that the ground state of QCD at finite quark density can 
behave as a color superconductor has attracted much attention recently |[l|-[16| . Of particular 
interest is the symmetry breaking pattern that may arise for Nf = 3 light quark flavors. 
As first shown by Alford, Rajagopal and Wilczek the diquark condensates can lock 
the color and flavor symmetry transformations (Color-Flavor-Locking or CFL for brief). In 
the limit of massless quarks, these condensates spontaneously break both color and chiral 
symmetries. The eight gluons become massive through the Higgs mechanism while eight 
plus two Goldstone bosons are leftover. The two octets of states are analogous to the 
octets of vector and light pseudoscalar mesons in vacuum and Schafer and Wilczek ^ have 
conjectured that there might exist some sort of continuity between the properties of QCD 
at zero density and in the CFL phase. The two other massless excitations are the Goldstone 
modes associated with the spontaneous breaking of the baryon number {U{1)b) and axial 
{U{1)a) symmetries. Because of the axial anomaly, the latter is not a true symmetry of 
QCD. However, because instantons effects are small at large densities, the associated mode 
can be treated as a true Goldstone mode as a first approximation. Although there are 
significant differences, to which we shall come back later, the situation at high density is 
analogous to considering the limit of large numbers of colors in vacuum |T^, in which 
U{1)a breaking effects are l/N^ suppressed. 

For energies which are small compared to the gap of the superconducting CFL phase, 
the dynamics of these Goldstone modes is most conveniently described with the help of an 
effective lagrangian. As suggested by the symmetry breaking pattern in the CFL phase at 
large densities, this effective theory is analogous to Chiral Perturbation Theory (xPT) in 
the large Nc limit of QCD in vacuum |]18|. For large densities, the leading order effective 
lagrangian invariant under U{3)l x U{3)r flavor symmetry has been constructed in some 
recent works [|T9| - |22| , |2^ , ^ . In particular. Son and Stephanov have shown how the pa- 
rameters of the lagrangian can be computed at large densities by matching to the underlying 
microscopic theory. In the present note we make a first attempt to extend these works to 
lower density regimes, taking into account the effects of instanton-induced interactions. By 
using the power counting rules of xPT we construct the effective lagrangians up to order 
E'^ in an energy expansion. We make a particular emphasis on the meson spectrum and 
compute their masses as function of the quark chemical potential /x. For moderate densities, 
II < 10^ MeV, the gauge coupling grows large and instanton effects become important. We 
give numerical values for the masses of the mesons simply assuming that the analytical ex- 
pressions for the gap and condensates, which have been computed at weak gauge coupling, 
also hold at strong coupling. Our calculations suggests that some of the Goldstone modes 
cannot exist as low energy excitations for values of /i < 1000 MeV, as their masses are very 
close to the unstability threshold 2 A, where A is the gap in the CFL phase. These conclu- 
sions depend on the value of the color superconducting gap as estimated in the literature, 
and of a number of approximations that we have made. 

In the next section, we write down the general low energy effective action for the pseu- 
doscalar Goldstone modes up to in a low energy expansion and estimate the couplings 
that are relevant for the meson mass spectrum. In Sec. ^ we give some numbers for the 
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meson masses at low densities, jj < 2500 MeV and finally draw some conclusions. 



II. EFFECTIVE CHIRAL LAGRANGIAN IN THE CFL PHASE 



We follow Gatto and Casualboni and Son and Stephanov to construct the 
effective lagrangian for the low energy excitations of the CLF phase of QCD. The ground 
state is characterized by the two diquark condensates 

where a, 6, c denote color indices, while refer to flavor ones. Under an S't/(3)c x 

SU{?>)l X SU{3)r the condensates transform as 

X -> UlXUI , Y ^ UrYUI , (2.2) 

and as 

X e^'" e^"^X , Y e~^'" e^'^Y , (2.3) 

under U{1)a and U{X)b transformations defined as 

^ e^^^+'^Vi , ^ c'^-^+^Hr . (2.4) 

One can factor out the norm of the condensates and consider the unitary matrices X and 
Y. The slow variations of the phases of these matrices then correspond to the low energy 
excitations. Altogether these are 9 + 9 = 18 degrees of freedom: 8 will be absorbed by the 
gluons through the Higgs mechanism, which leaves 10 true low energy excitations. At low 
energies the gluons decouple from the theory, as they are heavy degrees of freedom. It is 
convenient to collect all the Goldstone modes in the unitary matrix 

S = XY^ , (2.5) 

which is singlet of SU{3)c and U{1)b and transforms as 

S ^ e^^" UL^Uji , (2.6) 

under SU{3)l x SU{3)ji x U{1)a- This leaves apart the low energy excitation that emerges 
from the spontaneous breaking of baryon number symmetry U{1) b- Because baryon number 
is an exact global symmetry of QCD, this Goldstone mode is always massless in the CFL 
phase, independent of the quark masses. As our focus here is on the effect of chiral symmetry 
breaking by finite quark masses, to simplify our discussion we will simply drop this degree 
of freedom in the sequel. 

To proceed we parametrise the unitary matrix S as 

E = exp(tj] , (2.7) 
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where $ = ^^T"^ , A 
of SU{3) and = 



= 1, ... 9. The for A = 1, . . . , 8 are the Gell-Mann generators 
2/3 13, all normalized to Tr(T^T^) = 26^^. We use the same 



nomenclature as in vacuum for the nonet of Goldstone modes, 



1 

V3 



f Vo 

JVO 
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1 
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(2.8) 



In (|2.8| ), /,r and are the decay constants at finite density respectively of the octet of 
mesons and of the singlet 77°. A priori there is no reason for these to be equal. By matching 
to the microscopic theory, Son and Stephanov [^] have found 



J TT 



21-81og2 /i^ 



3 fi' 



f2 

18 27r2 ' 4 27r2 



(2.9) 



where n is the quark chemical potential. These expressions are valid at large densities [i.e. 



small gauge coupling). From (|2.9| ), the ratio fn/ frjo ~ 1-07 is close to one. This result is 
reminiscent of the OZI rule often invoked to hold in vacuum. At smaller densities, the 
gauge coupling grows large and the ratio could significantly depart from unity because of 
instanton effects. 

The quark mass term in the microscopic lagrangian 



AC = -^PlM^Pr + h.c. 



(2.10) 



breaks chiral symmetry. Its effects can be introduced in the effective lagrangian by treating 
the mass matrix as an external field with vacuum expectation value 



M = diag(m„, mrf,m^ 



(2.11) 



where m^, and refer to the up, down and strange quark masses, respectively. The 
spurion field then transforms as 



(2.12) 



under SU{3)l x SU{3)r x U{1)a. 

To take into account the effects of the U{1)a anomaly, we also allow for the presence of 
the 6 term in the microscopic lagrangian. 



(2.13) 



As for the quark mass term, we will treat 9 as an external field with vanishing expectation 
value which, to account for the variation of the quark measure, transforms as 



9^9- 2Nfa = 9 -6a 
under U{1)a transformations. Then any arbitrary function of the combination 



(2.14) 
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(2.15) 



is invariant under SU{3)l x SU{3)r x U{1)a transformations. 

With these ingredients, we are almost ready to construct the low energy effective la- 
grangian. One last issue is that of power counting. In QCD in vacuum, because ~ m^, 
the expansion is in powers of the pion external energy or momenta and quark masses and 

^ p2 ^ Similarly, at large Nc, U{1)a breaking effects are counted as E"^ ~ l/^c 
because M^^ ~ 1/A^c- In the CFL phase of QCD on the other hand, the power count- 
ing depends very much on the density. At finite densities, instantons effects are screened 
oc (A/yu)", where A ^ 200 MeV is the scale of QCD, fi is the quark chemical potential 
and a is some positive constant which depends on the process under consideration. This 
has two consequences at large densities. The first is that U{1)a is essentially a good sym- 
metry and the t]q is approximately massless in the chiral limit. Then, because the two 
condensates break chiral symmetry only through the intermediate of color transformations, 
there is an approximate Z2 x symmetry which acts independently on the left and right- 
handed quarks . At high densities, this implies that the leading contribution to the mass 
of the Goldstone modes is 0{m'^). The natural power counting rule at high densities is 
then E"^ ~ because ~ m^. At the moderate densities that could eventually be 
of interest for heavy ion collisions or for neutron stars, instantons effects are likely to be 
non-negligible. As x symmetry is broken to the diagonal Z2 by instantons, a meson 
mass term ~ ruq is allowed [0]. If instanton effects are dominant, the power counting is 
that of vacuum, E"^ ~ rUq. For convenience, we will adopt the counting rules of xPT as in 
vacuum and comment when necessary about the differences that arise in the CFL phase of 
QCD. 

A final remark concerning the region of applicability of the energy expansion. In vacuum, 
the low energy expansion is valid for E"^ < f"^. In the CFL phase however, ~ /i, which is 
independent and much larger than the gap (A ^ /i) at large densities (see (|2.9|)). However, 
the effective theory must break down for E ~ 2 A, which is the energy necessary to excite one 
quasiparticle pair out of the superconducting ground state, and the low energy expansion 
is only valid for the more restrictive range E"^ ~ < 4A^ f|. 

At order E"^, the most general lagrangian compatible with the symmetries of the con- 
densates is 

£2 = Vi{X)Ti (S^SS^St) + [V2{X)Ti (MT)) e'' + h.c] 

+ V?XX) (Tr Sa^St)' + V^{X) (TrSa^St) 3^6 + V^{X)d^ed^'e . (2.16) 



* There is a formal analogy between this behavior and that of xPT in vacuum in the large Nc limit. 
A priori the expansion is valid for E < Jt^. But because ~ NcA"^ grows large as A'c — > 00 and 
as something is bound to happen for ~ A ~ 200 MeV, which is the scale of quark confinement, 
the low energy expansion is limited to £^ < A. 
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We have written ( |2.16| ) using a compact notation. Since Lorentz invariance is broken at 
finite density (only a 0(3) symmetry is preserved) all the four- vectors should be split into 
temporal and spatial components, and the functions Vi multiplying those spatial or temporal 
components are not forced to be the same by symmetry arguments. For example, the first 
term in (12.161) should read 



V,{X)Tt (d^Ed^E^) ^ Vi,t{X)TT (aoEaoSt) - V^i,,(X)Tr {d,Ed0) • (2.17) 

The same situation occurs for the remaining terms and functions V3, V4, V5 of X. All the 
couplings in ( ^.l(j| ) can be a priori arbitrary functions of X ( ^.15| ). At large densities, U{1)a 



symmetry breaking effects are exponentially suppressed and the couplings only depend on 
the chemical potential /i. 

To reproduce the standard normalization of the meson kinetic terms, we impose 

Vi,m = f . (2.18) 

The ratio Vi_s(0)/Vi^t(0) = is the velocity squared of the Goldstone bosons. At large 
densities. Son and Stephanov |2^ have found that v is equal to the speed of sound l/Vs 
for all the low energy modes, including the baryon Goldstone mode. 

The operators in (|2.16|) are not all independent. The last three terms can be transformed 



into each other with a field redefinition, '^0//% ~^ Vo/ frjo + f^d- Using this freedom, we can 
choose to set V4(0) to zero. With this choice, the last operator becomes irrelevant for the 
meson spectrum and can be discarded. 

The operator that is left, with coupling V^^X), contributes to the difference between 
and /^Q. At high densities, using 
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^3,(0) = ^^j^ ^ OMfl (2.19) 

which is small compared to At moderate densities, V3 1 could receive large contributions 
from instantons as is manifest from the mixing with V4 and V5. 

Consider now the mass term in ( |2.16|) . This term is analogous to the leading mass 
term in xPT in vacuum. The only difference is the occurrence of the phase 6, which 
is absent in vacuum. This is because at zero density the condensate that breaks chiral 
symmetry is {4'l4'r) which transforms like under SU{3)l x SU{3)r x U{1)a- In ( p.l6| ) 



the presence of the 6 in the effective lagrangian is the trademark of a one instanton effect; 
in the CFL phase, a one instanton process can be saturated by closing its six external quark 
legs with the insertion of one left-handed diquark, one right-handed diquark and one chiral 
condensate, ~ (V'lV'l)(^-R''/'r)(V'-R^l)- As in vacuum, a non-zero chiral condensate (ipip) 
leads to oc mq{il)il)). Instanton effects are small and can be reliably computed at high 



densities /i ^ A. In particular, Schafer |]I5[ has obtained 

{ifjip) = {uu + dd + ss) 

" ''-G(/x)0i(/i) , (2.20) 



fi^ \ 3V27r 18^, ^ ,2 



2n^J g{fi) 5 
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where g{fi) is the gauge couphng estimated at the scale fi using the one- loop beta function. 
The factor G{n) is the one instanton weight integrated over all instanton sizes p, which are 
peaked around p ~ at finite density, 



G{p) ^ 0.26 A-^ ((/?o log(/i/A))'^^ f - j . (2.21) 

We have used the running of g{n) at one-loop, A = 200 MeV and (3q = 11/3 Nc-2/3 Nf = 9 
is the first coefficient of the QCD beta function. Finally, (/'^(/i) ~ ~ {i'R'ipR) is the 

diquark condensate in the CFL phase of QCD, 

M^^) - 2 f ^1 A(p) , (2.22) 



2vrV V ^?(/i) 



and 11,5,13 



A(^)=6;512i-'(2/iV^)=/'^g(^)-=^expf--|^l , (2.23) 

is the color superconducting gap. The factor 6q is unknown but expected to be C(l) 
TB| , |2B|P7| . At large densities /i ^ A, instantons are suppressed, (iJip) goes to zero and 



the mass term of ( p.l6| ) vanishes. At the moderate densities that could be of interest for 
heavy ion collisions or neutron stars, this instanton effect is however not negligible. Using 
the Gell-Mann-Oakes-Renner relation, and (ipip) = {uu + dd + ss) = 3{uu), we find 

V2{0) = -U^P^P) . (2.24) 

At order E^, there are a few more operators^, 
£4 = Ki{X) [Tt (^d^J:d>'E^)f + K2{X)Tt (S^SS^St) Tr (d^'Wi:^) (2.25) 

+ {K4{X)Tt (S^SS^St) Tr (AlEt) ^ | ^ ^K5{X)Tt (a^Ea^S^A^S^) e^^ + h.c.} 
+ {Kq{X) det(S) Tr (tWE^) Tr (A^S^) + h.c.} + K7{X)Tt (A^S^) Tr (A1"^s) 
+ {Ks{X) det(S) Tr (MS+A^St) + h.c.} . 

Again, we have used a compact notation for the terms in i^^j, i = 1 . . . 5 which should be 
split into temporal and spatial components. 

Like the term V2 in (|2.16|) , the terms K4 and i^s depend explicitly on 6 and so are 
exponentially suppressed at high densities. Similarly, at large chemical potential n, the 



§At this order, one could also consider adding a Wess-Zumino-Witten term |1T 
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remaining functions Ki{X) reduce to constants Ki{X = 0; /i). In this limit, the mass pattern 
of the Goldstone bosons will be determined by the couplings Kq, Kj and Kg in (|2.25|) . As 
first shown by Son and Stephanov at large densities, these coupling constants can be 
computed by matching to the underlying microscopic theory. 

In principle, a systematic and non-ambiguous strategy to compute the coefficients of the 
effective theory is to use the background field technique p8[, introducing external sources 



and symmetry breaking order parameters and integrating out the quark fields. At high 
densities, gluon exchange is suppressed and this amounts to a quark one-loop calculation, 
which would fix the coefficients of all the operators to arbitrary order in the meson fields. 
A slightly different strategy has been advocated in |]22|: set the meson fields to zero, S = I3 



in the operators of (|2.25| ), and compare the shift in ground state energy induced by non- 



zero quark masses in both the effective and microscopic theories. This approach is a priori 
ambiguous, as different operators could contribute to the shift in ground state energy. In 
practice, to order E^, enough constraints can be derived to completely fix the couplings Kq, 
Kj and Kg. The diagrams in the full microscopic theory are those of Fig.l. The rhs diagram 
can only contribute to Kq and Kg. The Ihs diagram can contribute to Ky. These couplings 
can be fixed by considering two different quark mass patterns, 

Mi = ml^ and A<2 = diag(0, 0, m,) = ^ + ^ . (2.26) 

In the effective theory, these respective choices lead to the following shifts in ground state 
energy density 

Aei = - {(WfsTg + h.c.) + WiTr + (Sm^/sTg + h.c.)] , (2.27) 
A£2 = - [{jtiIKq + h.c.) + mlKj + {mlKg + h.c.)] . (2.28) 

A special feature of a quark mass term is that it couples fermions and antifermions. At 
finite density, quark mass effects are then suppressed by the necessity to excite antiparticles 
with characteristic momentum ~ 2/i p^ , ^ . A simple way to see this is to consider the 
effective theory for the fermion excitations near the surface of the Fermi sea. In momentum 
space, the free fermion lagrangian is 

C = ^{Yp^ + fi-f^)^ - rmp^ . (2.29) 

If we consider the fermion mass term as a perturbation and decompose the fermions field i/j 
into positive and negative energy components and ip^, using the projectors 

A± = ^(l±a-q) , (2.30) 

and A^ip± = il)±, the lagrangian becomes 

^ = i^Uqo -\q\+ /i)V^+ + i^U<io + \q\+ f^)i^~ - m(V'U- + • (2-3i) 

This shows that the fermion mass term couples particles and antiparticles. If we integrate 
out the antifermions, we get at leading order 
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2 

m 



£^4(9o-g|i)^+ (2.32) 

where ~ |^ — This shows that quark mass effects are suppressed at high densities, as 
one could have expected. Note that the modified mass term oc m^/ /i has the same structure 
as a chemical potential term. This is a relevant operator which modifies the shape of the 
Fermi surface p9 |. 



The calculation of the l.h.s. diagram of Fig. 1 gives 

^ = 3^ . (2.33) 

to 0{ ^4 ), where A stands for generic gap and/or antigaps. This diagram simply gives a 
trivial shift of the vacuum energy, which in the effective theory corresponds to a constant 
(independent of S) invariant operator 

AC oc Tt{M^M) , (2.34) 

and, consequently, Kj = 0. 

Estimates of the r.h.s. diagram of Fig. 1 indicate that [ p^ 



AA 

i^6,8~^log(A//i) , (2.35) 
A* 

where A is a generic expression for a quark gap and A is an antigap. The dependence on 
antigaps arises from the structure of the r.h.s. diagram with both particles and antiparti- 
cles propagators. Unfortunately, the expression of the antigap is not known. Preliminary 
estimates have revealed that it is gauge-dependent |]12|. It is however a physical quantity 
(pole of the antiquark quasi-particles and holes), which should be gauge-invariant on quasi- 
particles mass-shell. At weak coupling, the antigap is presumably much smaller than the 
gap. Here, we will assume that Kq^s ~ at large and moderate densities. 

To be complete we should also include the contribution of instantons to the mass of the 
singlet meson tjq, which enters the effective lagrangian through a function of X alone, 

Co = -Vo{X) , (2.36) 

This piece is as it involves no derivatives of the meson fields. Expanding to second order 
in X gives 



Mi = l^. (2.37) 



The constant Vq" (and more generally the whole function V{X)) could in principle be com- 
puted at large densities using instanton calculus. Unfortunately, the result is not known. 
However, at large to moderate densities we might expect this contribution to be suppressed. 
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For three flavors, the two remaining legs of a one instanton contribution to M^^ can only 
be closed with the insertion of either a quark mass term or a chiral condensate {^'ip)- 
the chiral limit, the latter only arises through another instanton process. On dimensional 
grounds we would expect 

Ml^^i'G{Mn)\ (2.38) 

in the chiral limit, but from such a rough estimate, it is not reasonable to infer whether M^^^ 
ever grows large at low densities. 



III. MESON MASS PATTERN 



In this section we give numerical estimates of the pseudoscalar mesons as function of the 
quark chemical potential. 

At very large densities, instanton effects are suppressed: meson masses are linear in the 
quark masses and presumably very small. At lower densities, a non-zero {ipip) condensate 
induced by instantons introduce contributions to the meson masses which are proportional 
to m^^. Here, we will estimate this last effect as it is dominant in the low density regime 
of the theory. We will use the expression of the chiral condensate and the gap as computed 
at weak coupling ( p^.20D and we will assume that the corrections are small even at moderate 
densities /i ~ 600 MeV. 

The masses of the charged pions and kaons deduced from (|2.16|) read 



-pr- {rriu + md) 



2V2 

Mj^± = — - {rriu + rus) 



J n 



(3.1) 



In the limit of no U{1)a breaking effects, the neutral mesons 7r'^,?78 and r/o are strongly 
mixed. Their mass matrix reads 



/ 



2V2 

J TT 



niu + md 



2(m.n-mrf) 

fir V2{niu~ma) 
V fvo 



JjL. 

fvo 



V3 
3 

\/2{mu+ma-2ms) 
V3 



U V2{mu-md) ^ 
/'70 \/3 
/tt V2(mu +ma -2ms) 
fvo V3 
fl 2(rnu+md+ms) 



(3.2) 



■"70 



We have assumed that the contribution from the quark masses and the chiral condensate 
is parametrically larger than the two-instantons contribution to the mass of 779. As the 
latter effect, if important, would increase the mass of r^o, our numerical estimates should be 
considered as lower bounds on the mass of rj' . Neglecting the difference between /^r and /^q 
we get the mass eigenvalues characteristic of ideal mixing. 



Ml 



2Vo 



f2 
J n 



2 2m„ , 



Mid 



f2 

J TT 



2m„ 



Ml 



f2 

J TT 



■ 2mo 



(3.3) 



To give an idea of the numerical values of the meson masses at /i ~ 600 MeV, we take 
mu = 4 MeV, = 7 MeV and = 150 and evaluate the value of the chiral condensate 
for A = 200 MeV and 6' = 1, to find 
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M^i ^ 30 MeV , M|± ^ 112 MeV , M|,o,^o ^ 113 MeV , (3.4) 
Ml^ ^ 26 MeV, ^ 34 MeV , Ml ^ 71 MeV . (3.5) 

A more accurate analysis of the meson masses would require to take into account terms 
which are proportional to m^. However, in the low density regime, and due to the smallness 
of Kq and Kg, these represent negligeable corrections to the above estimates. 

Because the expression for the chiral condensate involves the square of the gap, which 
is only known up to a factor 6q assumed to be 0{1) |lTT|,^, p^ , |26| , |27| , there is some theoret- 



ical uncertainty in the meson masses at low densities and the values quoted are at most 
indicative. At this order, this uncertainty can be cancelled by considering the mass over gap 
ratios. These ratios are independent of the gap and, moreover, tell us whether the mesons 
can exist as low energy excitations in the CFL phase at chemical potential fi. The ratios 
of the charged and neutral meson masses over two times the gap are shown in Figs. 4 and 
5. For the strange particles, the ratio is larger than one or very close to the instability 
threshold unless /x > 1 GeV, which suggests that these particles do not exist as stable low 
energy excitations at low densities. 

IV. CONCLUSIONS 

We have considered the effective lagrangian for the Goldstone modes in the Color-Flavor- 
Locking phase of QCD at high densities, up to in a low energy expansion E < 2A. A 
tentative analysis of the meson mass pattern that emerges from this lagrangian, including 
instanton effects, suggests that the kaons and the pure strange neutral meson may be absent 
from the spectrum for moderate densities. Their masses are significantly larger than two 
times the gap for /i < 1 GeV. For the pions, the situation is much better as there their 
masses are significantly smaller than the instability threshold for all densities. 

Obviously, there is room for improvements. Our approach has been far from systematic 
and we have made some ad hoc assumptions. Although the effects are likely to be small at 
moderate densities, it is of much interest to gain a better understanding of the quadratic 
quark mass effects plJSH] . A calculation of the contribution of instantons to the mass of r^o in 



the chiral limit would also prove to be most useful. If a spontaneously broken, approximate 
U{2)l X U{2)r flavor symmetry survives in the CFL phase of QCD at moderate densities, it 
could also be interesting to determine the potential for the r^o, Vo(X). The effective theory 
for the low energy excitations should be very much analogous to the large N^. lagrangian of 



Di Vecchia and Veneziano 31 and Witten 32 
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FIG. 1. Quark loops in the microscopic theory. M (M^) stands for one insertion of the (con- 
jugate) quark mass matrix. The Ihs diagram is Hke in vacuum. The rhs diagram arises only in 
presence of a diquark condensate. 
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FIG. 2. Color gap (in MeV) as a function of the quark chemical potential fi (in MeV) for I'q = 1. 
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FIG. 3. Plot of (in MeV) as function of the quark chemical potential fi (in MeV), 

where {'tpi/j) = {uu + dd + ss) is the chiral quark condensate induced by instantons. 
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FIG. 4. Ratios of charged pions (lower curve) and kaons (two upper curves) masses over two 
times the gap as function of the quark chemical potential (in MeV) . 
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FIG. 5. Ratios of neutral meson masses (Mu„ < Mj^ < Mgg) over two times the gap as functions 
of the quark chemical potential (in MeV). 
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